Resonant formation of strongly correlated paired states in rotating Bose gases 
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We propose increasing the fractional quantum Hall gap of a rapidly rotating Bose gas by increasing 
the interatomic interactions via a Feshbach resonance. The generation of molecules by the resonance 
causes pair correlations to grow throughout the system effecting the ground state. By an extension 
of the usual Chern-Simons theory, built of composite atoms and molecules, we are able to account for 
these resonance effects. We find that the resulting ground state evolves from a Laughlin wavefunction 
to a unique paired wavefunction as one approaches the resonance. 

PACS numbers: 03.75.-b,74.43.-f,71.27.+a 
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Recently, there has been considerable interest from 
both the experimental and theoretical perspective in the 
behavior of rotating Bosegases confined to an effective 
two-dimensional space 0,0111 El 01 El- One goal has been 
to understand and to create strongly correlated states, 
such as the Laughlin state generated by the fractional 
quantum Hall effect (FQHE), within ultracold atomic 
gases. A major challenge to realizing such a state is 
the need to reach extremely low temperatures in order 
to resolve the lowest Laughlin state from the first excited 
state. The size of this energy gap, however, is directly 
related to the strength of interatomic interactions, so by 
increasing the interactions it would seem natural that one 
could increase the gap, making the system more accessi- 
ble to experiment. 

Atomic systems can now be created which allow the 
microscopic interactions to be dynamically tuned 0, @ • 
Feshbach resonances have proven very successful at doing 
this, allowing one to tune the interactions by adjusting 
the resonant detuning v (for example, by varying a mag- 
netic field), and would seem an excellent tool for increas- 
ing the gap within a FQHE system. To account for the 
full effects of the resonance, however, we cannot simply 
scale the mean-field energy but must incorporate the en- 
tire resonant structure into our model. This means that 
we must include the process of molecular formation gen- 
erated by the introduction of a bound state within the 
open channel of scattering states 0, 0, 0] . 

By introducing a bound state, however, we not only 
modify the relative interaction strength, but we also in- 
troduce a physical mechanism for generating pair corre- 
lations between particles (see Fig. In the context of 
two-dimensional condensed matter systems, such a mech- 
anism, although arising from a very different source |l2| , 
can have a significant effect on the ground state proper- 
tics. Before we can study the resonant behavior of the 
gap, we must first understand the effect of resonant inter- 
actions on the ground state wavefunction. We will find 
that as we approach the Feshbach resonance the Laughlin 
state transforms into a unique, strongly correlated state. 

We begin by writing down an effective Hamilto- 
nian, in second quantized form, for a resonant gas of 
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FIG. f: The Feshbach resonance pairing mechanism is il- 
lustrated by the above Born-Oppenheimer curves. Pairs of 
composite atoms composed of single atoms and an associated 
number of quanta of angular momentum (represented by the 
arrows) approach each other within an open channel poten- 
tial of background value U{x, x'). They may form a compos- 
ite molecule due to the presence of a closed channel bound 
state, at a detuning v from the scattering continuum, which 
is coupled to the open channel with strength g(x, x'). 



Bosons of mass m rotating in two dimensions with stir- 
ring frequency £1 approaching the trapping frequency 
uj, i.e. u; — O — > + : 



H= /d 2 x^(x) 



-(V-.AW)* 



(1) 



i \d 2 x'\d 2 x ^t (x )^t (x ' )r/(X)X / ) ^ a(x / ) ^ (x) 



d 2 x $„(x) 



d 2 x' ld 2 x 



4m 



(V - 2iA(x)) 2 



^(^ L )5(x,x')^ a (x)^ a (x') +H.c. 



Here m (x) , ij) a m (x) are the creation and destruction 
operators for atoms and molecules which satisfy the 
commutation relations [V' 1 (x), ^|(x')] = (x,x')5i^, 
where 1 , 2 £ {a, m}. We define the two-dimensional vec- 
tor potential A(x) = (muy, —mux), C/(x, x') is the two- 
particle background scattering potential, <?(x, x') is the 
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resonant coupling between the open and closed channel 
potentials, and v is the detuning of the open channel con- 
tinuum from the level of the bound state in the closed 
channel. 

Recently, it has been noted that there is a direct 
mapping between the Hamiltonian for a rotating two- 
dimensional gas, the first two terms of Eq. (JTJ , and that 
of the FQHE |4j . The analog of the magnetic field is real- 
ized by the angular rotation and the Coulomb interaction 
is replaced by the two-particle scattering. Therefore, if 
we were to neglect resonant effects we would expect the 
many-body ground state to be described by the Laughlin 
wavefunction ll.'l 



tf L (xi,x 2 ,...x JV ) = Y[(z l -z ] ) 2 Y[exp(-\z k \ 2 /2) , (2) 

i<j k 

where the products run over the indices k — 
(1, 2, ...N) at position z = x + iy for TV particles. We will 
see that the inclusion of the resonant terms in Eq. 
can significantly modify the form of Eq. (|2J due to the 
growth of two-particle correlations. 

We approach this problem by an extension of the 
Chern-Simons theory [14j which allows us to develop a 
mean-field theory for the rotating system that has re- 
moved the complications of the associated rotation. This 
is done by constructing a composite particle composed 
of the original particle and an artificially attached num- 
ber of flux quanta. The composite particle is designed so 
that the attached flux quanta cancel the total rotation 
of the original system leaving a system of non-rotating, 
interacting composite particles. For the resonant system, 
the composite particles can formally be obtained by the 
following transformation: 



0a : m(x) = CXp 



- q j d 2 x'6{x - x') 
x [p a (x') + 2p m (x')] 



^a,m( x )> (3) 



where p a {x') and p rn (x') are the atomic and molecular 
spatial densities, respectively, and 6(x — x') is the topo- 
logical phase. To guarantee that the resulting composite 
particle is a Boson, with no loss of generality we set q = 2. 

In the composite picture, we replace the tp operators 
in Eq. with the corresponding composite operators 
tp and introduce the statistical Chern-Simons field a(x) 
through the gauge transformation 



A(x) -► A(x) + a(x) 



(4) 



These modifications generate the Hamiltonian formula- 
tion of our composite atom/molecule theory. The com- 
posite picture can be shown to be equivalent to the single 



particle picture of Eq. QJ . 

We will now shift to a functional representation of the 
composite atom/molecule system to clarify the resonant 
modifications of the Chern-Simons theory and then re- 
turn to the Hamiltonian formulation to derive the ground 
state wavefunction. Defining the action within one tem- 
poral and two spatial dimensions 



S = Jd 3 x VlWiBovA*) ~ [ dx oH, (5) 
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we generate a Chern-Simons term which couples to the 
statistical vector field a M (x) 



d 3 x— e^ x a^x)d u a x (x), (6) 



where the indices of p, v, and A run over the three di- 
mensions (0, 1, 2) and the summation convention over re- 
peated indices is invoked. We have also introduced the 
antisymmetric tensor e 012 = 1. To simplify the follow- 
ing calculations we assume contact interactions of the 
form U(x, x') = US(x, x') and g(x, x') = gS(x, x'). Any 
complications arising from this replacement of the true 
potentials with contact potentials should be remedied as 
explained in Ref. 01 • We next perform the lowest order 
variation of the action. 

Varying with respect to the zeroth component of the 
gauge field, dS/da — 0, reproduces the Chern-Simons 
condition 



Vxa(x)|^ = -47r(|^ a (x)| 2 + 2|^ m (x)| 2 ) 



(7) 



Equation (JJJ is a statement of Gauss's law for the statis- 
tical gauge field associating an even number of rotational 
flux quanta with each particle. This relation is simply a 
restatement of our choice of quasiparticle. 

Since we will be interested in the ground state proper- 
ties of the atom/molecule system, let us assume that the 
fields f a (x) and tp m (x) are uniform. By minimizing the 
action with respect to the atomic and molecular fields, 
i.e. dS/difa = 0, dS/dip m = 0, we generate the following 
constraint equation for the molecules: 



<Pt, 
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2[u + | A + a| 2 /m] 



(8) 



Equation © allows us to eliminate the molecular field 
from the theory and arrive at a self-consistent relation- 
ship for the gauge field 



|A + ar = \ U 



4(V + |A + a| 2 /m) 



2m\p a \ 2 . (9) 
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Equation |j5Jl is the usual result relating the gauge field 
to the background density only now it is dependent upon 
the detuning from the resonance. 

We now switch back to the Hamiltonian form of our 
theory to derive the ground state wavefunction. After 
Fourier transforming the composite form of Eq. Q by 
substitution of the field operators <^ a (x) = Ek"k f 



ik-x 



and y? m (x) 



Ek^kE 



;k-x 



we follow the usual 
Hartree-Fock-Bogolubov (HFB) approach to construct a 
quadratic Hamiltonian which accounts for the lowest or- 
der pairing. As before, we assume contact interactions 
and now make the additional assumption that we may 
neglect the excited modes of the molecular field keeping 
only the lowest condensed mode 6k = &o- The resulting 
Hamiltonian for the composite system can be written in 
the form 



(10) 



k#0 



H° is composed of all terms of less then quadratic or- 
der in the operator d k , we define a column vector A k — 
(d k , a_ k ), and M k is the self-energy matrix. For our pur- 
poses we need only concern ourselves with the structure 
of the second term in Eq. I)l()|l. Here the self-energy ma- 
trix is expressed as 



where Ti. contains the ground state contribution to the 
energy and the excitations are given by the spectrum of 
frequencies Wk = \/E^ — |A| 2 . 

Since there are no quasiparticles present in the ground 
state \gs) , which is what one would expect from an 
interacting bosonic system at T = 0, the ground state 
must satisfy the condition Skiffs) = 0. Substitution of 
Eqs. (|14f> and i|15|) for the quasiparticle operators result 
in the relation 



(E k + uj k )a k \gs) = —Aa!_ k \gs). 



(17) 



Because d k and d k are canonically conjugate variables 
there is no loss of generality in making the replacement 
d k — > <9/i9d k [lfij . This converts Eq. (|17|l into a simple 
differential equation for the ground state with the solu- 
tion 



\gs) = exp 



L k 



-A 



- w k 



« f -k4 



|0>. 



(18) 



To derive the many-body wavefunction we must now 
move from second to first quantization. The relationship 
which links the second quantized ground state \gs) with 
the first quantized wavefunction $cb can be written for 
an even number of noncondensed particles 2N as 



Mv 



E k A 

A* £_ k 



*Cb(xi,X 2 , ...X 2 jv) = (0|(^ a (x27v)---^a(x 2 )^a(xi)|ffs), 

(19) 



with the additional definitions for the diagonal and off- 
diagonal terms 



EZ + U(\<p a \ 2 + n) 



A = U(cpl+p) + g<p r , 



(12) 
(13) 



Equations (|12fl and l|13|) are expressed in terms of the 
pairing-field p — ^ k , (d k ,d_ k ,), the normal-field n — 

Ek' (^L^k')' ano - ^k ^ s the effective kinetic term which 
contains the contribution from the gauge field ^4(x). 
Eq. (|10l) can be rewritten in terms of the quasiparticles 



where it should be noted that $cs is the full-many body 
wavefunction for the composite Bose particles. If we are 
able to assume that E k + w k 3> A , an assumption which 
will remain valid as long as we are not too close to res- 
onance, we may truncate the power expansion of the ex- 
ponent in Eq. I|18|l . This results in the composite boson 
wavefunction 



*Cs(xi,X 2 , ...X 2 a0 = SW>12^34— i>{2N-l)2N) (20) 

comprised of a symmetrized product S of paired wave- 
functions 



a - ^ (Ek + ^ + |Ap ( (Ek + - K + A °'-) (14) (21) 

(^(E k + w k )a_ k + A t a k ^) (15) If we were dealing with a system of fermions, equa- 
tion (|2*U|l would be antisymmetrized and would result in 



V(^k + ^k) 2 + |A 



symm 

which result in the diagonal Hamiltonian a Pfaffian wavefunction [Tfl. Here, because of the statis- 

tics of the particles, we generate a bosonic analogue to 
this result. The many-body wavefunction for the bare 
Ti = Ti. + Wk« k a k , (16) particles can now be extracted from the composite wave- 

k ^o function [l8| resulting in 
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= *cs(xi,x 2 ...,x 2 Ar) X * L (x x , X 2 . .. , X 2 N ) (22) 

which is a product of the composite particle wavefunction 
of Eq. (|2U[1 and the Laughlin wavefunction of Eq. @ . 

Equation (|22[1 is the final result for the ground state 
wavefunction of the resonant rotating Bose system. This 
result has important consequences for the generation of 
the FQHE within a resonant atomic gas. It would imply 
that an increase of interparticle interactions by a reso- 
nant tuning of the interactions simultaneously generates 
an increase in 2-particle correlations between compos- 
ite particles resulting in a modification to the ground 
state wavefunction. As is clear from the form of Eq. (|20H . 
for large detuning from the resonance, corresponding to 
small pairing and molecular field, the many-body wave- 
function reduces to the Laughlin wavefunction. As one 
moves nearer to the resonance, however, the off-diagonal 
part of the self-energy matrix, A, grows. This results 
in an increasing modification of the many-body wave- 
function from the Laughlin wavefunction. The ability to 
tune a Feshbach resonance therefore allows for the direct 
study of this crossover from a Laughlin wavefunction to 
a paired wavefunction. 

A similar, yet distinct, paired wavefunction as in 
Eq. (|22) was found for electronic FQHE systems 0, 03 
This has been used to explain the previously un- 
resolved even denominator filling fractions which result 
from a pairing instability, such as the observed incom- 
pressibility of the 5/2 filling. In this straightfor- 
ward generalization of the Laughlin wavefunction would 
result in a symmetric wavefunction, violating the asym- 
metry of the fermions. However, the generation of an 
antisymmetric Pfaffian wavefunction which multiplies the 
generalized Laughlin state allows the overall ground state 
to be correctly antisymmetrized. For the bosonic system 
we have treated, the overall wavefunction must remain 
symmetric, so the corresponding paired wavefunction is 
symmetric in comparison to the antisymmetric Pfaffian 
wavefunction. 

In conclusion, because of the extreme diluteness of 
trapped atomic gases, if strongly correlated effects such 
as the FQHE are to be observed in these experiments, the 
interatomic interactions will most likely need to be reso- 
nantly enhanced. Feshbach resonances are a convenient 
method for increasing the interactions, but the molecular 
processes involved in such a resonant system force one to 
account for the effects of atomic pairing. The growth of 
pair correlations among the composite particles leads to 
a modification of the expected ground state for a rapidly 
rotating Bose gas. The new ground state wavefunction, 
which is generated by the Feshbach resonance, exists as 
a strongly correlated state unique to trapped Bose gases. 
The tunability of the resonance opens the possibility for 
the direct study of the crossover transition between the 



paired state and the Laughlin state. 

These results have important implications for the pro- 
duction of the FQHE within atomic gases since the use 
of a Feshbach resonance results in a state quite differ- 
ent from a rapidly rotating gas where the interactions 
are quantified by only a large scattering length. For in- 
stance, many of the observeable properties of the gas may 
be modified such as the density profile for both atoms 
and molecules and the nature of collective excitations. It 
should also be noted that the crossover transition we have 
discussed is only a part of a much more general crossover 
theory made accessible by the tunability of a Feshbach 
resonance. Although the methods presented here are in- 
valid close to the resonance, one could imagine extending 
these ideas to describe the resonant system as one passes 
from a gas of interacting rotating atoms, through the res- 
onance, to a system of tightly bound, rotating molecules. 
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